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Introduction 


Tc  connp^r’.rs  po 
phase  shift  may  be 
variable  type  phase 


\:er  and  increase  the  accuracy  to  nnioh  a  prescribed 
set,  digital  rather  than  the  usual  continuously 
shifters  have  been  suggested  for  use  in  phased  micro- 


wave  antenna  arriysa  Various  schemes  have  been  suggested  including  the 
use  of  magnetically  saturated  toroids  in  rectangular  t^vegvides  in  which 
the  phase  can  bo  shifted  digitally  by  the  reversal  of  the  magnetization  in 
the  toroid  thi'ough  the  application  of  a  current  oulse,  tfeny  cf  these 


arrangements  suffer  from  the  disadvantage  that  not  all  of  the  ferrite  is 
effective  in  producing  diffei'ontial  phase  shift  since  only  a  portion  is  in 
the  region  of  circular  polarization* 

Ferrite  toroids  in  round  waveguides  utilizing  the  circular  TE  mode 
seem  to  be  the  natural  mici’owave  structure  for  such  non-reciprocal  phase 
shifters  for  several  reasons^^\  First,  it  is  possible  to  place  all  of 
the  ferrite  in  a  region  cf  circular  polarization  so  that  the  entire  vol¬ 
ume  of  material  is  then  effective  in  producing  non-reciprocal  phase  shift. 
Secondly  it  is  possible  to  make  narrow  transverse  cuts  through  such  a 
structure  without  introducing  significant  reflection  so  that  a  number  of 


these  section!^  may  be  cascaded  to  allow  d,  c.  insulation  bet^/een  them;. 
Thirdly,  such  transverse  cuts  may  be  used  for  the  introduction  of  wires 


to  carry  the  svjitching  current  to  the  toroids.  Finally.,  a  thin  wire  may 
be  introduced  along  the  waveguide  axis  vdthout  aopreciably  changing  the 
fields.  Tl'.e  ferrite  magnetization  could  then  be  sidtched  by  a  current 
flovring  along  this  wire.  Transverse  cuts  in  this  center  conductor  can 


The  use  of  the  TE^^  mode  in  conjunction  wdth  circumferentially  rvuagnet- 

ized  ferrite  tubes  inserted  in  circular  cylindrical  guides  appears  to 
have  been  suggested  first  by  A.  G.  Fox,  S.  E.  Fuller  and  M.  T.  V7eiss, 
"Behavior  and  Application  of  Ferrites  in  the  I-Iicrowavo  Region",  Bell 
Sys.  Tech.  J.?  Vol .  PP«  dan.  1955”  The  use  of  such  struc¬ 

tures  as  digital  phase  shifter^,  was  proposed  the  second  author  v.hile 
he  v;as  at  Ldncoln  Laboratories,  see*  Solid  State  Resoarchj  Lincoln 
Laboratories,  ifess.  Inst,  of  Techn- ,  No.  1,  Sectior  T/  Q,  pp.  38,  I963. 


al30  be  niade  ’.cithout  disturbing  the  mode. 

The  usual  microwave  ferrite  materials  nay  not  be  suitable  for  such 
application  since  they  \d.ll  require  desirable  s^idtching,  as  well  as  good 
microwave  properties.  Such  materials  are  under  investigation  at  this  as 
well  as  other  laboratories, 

A  theoretical  investigation  is  presented  here  of  the  mode  structure 
and  the  differential  phase  shift  per  \init  length  for  the  cii'cular  cylin¬ 
drical  guide  containing  a  coaxial  ferrite  tube,  .Jhile  it  is  recognized 
that  in  practice  the  ferrite  toroid  will  be  of  finite  length  and  the 
ferrite  vdll  not  necessarily  be  completely  magnetically  saturated  it  is 
expected  that  this  study  vd.ll  yield  useful  design  data  such  as  the  dif¬ 
ferential  phase  shift,  optimum  toroid  placement  and  dimensions. 
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Theory 


We  will  first  consider  wave  propagation  in  a  circular  cylindrical 
guide  comnlet'^ly  filled  with  ferrite  magnetized  in  the  circumferential 
or  0  -  direction.  The  ccin"^lGtely  filled  gjide  is  treated  only  to  study 
the  charactei’  of  the  eigenfunctions  and  their  eigenvalues  and  it  is  not 
expected  (due  to  the  difficulty  of  synthesizing  a  uniform  circumferential 
magnetizing  field)  that  such  an  arrangement  will  have  practical  utility. 
The  analysis  wi.ll  then  be  extended  to  more  complex  structures  such  as 
tubes  or  rods  of  ferrite  placed  coaxially  within  circular  cylindrical 
guides. 

A  ferrite  magnetized  in  the  direction  of  the  0  -  coordinate  nay  be 
characterized  by  the  permeability  tensor 


where 


0 


1+X  0  jx 

0  10 
-jx  0  1+X 


TM  rri  ,  ^  ^  coYH 


(1) 


(2) 


and  Y  is  the  gyromagnetic  ratio.,  M  the  magnetization  and  H  the 
magnetizing  field.  V/e  will  restrl-'t  our  attention  to  those  TE  modes  that 
exhibit  rotational  symmetry.  Ass\iming  time  and  axial  variations  of  the 
form  and  e**^*^^,  respectively,  we  find  that  the  only  non- vanishing 

component  of  the  electric  field  must  satisfy  the  differential  equation 

or  r  ^ 


(3) 


where  r  is  the  radial  coordinate  and 


a»  -  px/(l  +  X) 

C  A«  -  (I4) 

o 

A*  =  [  (1  X)^’  -  x^3/(l  +  X) 

The  components  of  the  magnetic  field  are  given  by 


H 


+ 


X  1  d 
r  dr 


u 


(1-^X)  1  ^ 
2  r  dr 
Yo 


V 


where 


[(l+X)^-x^3  •  A  A 

0  0  o 


We  seek  solutions  to  (3)  subject  to  the  boundary  condition 


(5) 


(6) 


-  0,  r  -  b  (7) 

where  b  is  the  inner  radius  of  the  guide. 

Equation  (3)  is  recognized  as  Tricomi’s  form  of  the  Confluent  Hyper- 

(2) 

geometric  equation  ,  V/e  may  thus  constmict  solutions  using  Humbert 's 
or  Whittaker’s  notation  or  other  forms  of  the  Confluent  Hype rgeome trie 
functions^^^ ,  Since  the  tabulation  of  th''se  functions  and  their  eigen¬ 
values  is  not  complete  it  was  decided  to  construct  solutions  that  relate 
directly  to  Bessol  and  Neumann  functions.  This  was  prompted  by  the  fact 
that  as  o’  goes  to  zero  (3)  reduces  to  Bessel’s  equation.  Further,  as  a' 


^"■^Bateman  Con^^endia,  McGraw-Hill  1953>  Vol.  1,  p«  2$]., 

(y) 

''Morse,  P.  M,  and  Feshbach,  H.,  Methods  of  Theoretical  Physics, 
McGraw-Hill  19?3,  Vol,  1,  pps,  60li-6l9, 


xS  proportional  to  tlis  product  Px  have  tuat  the  eigenvalues  must  reduce 
to  those  of  the  Bessel  functions,  both  at  cut-off  P  0,  and  as  x  goes 
to  Of  i.e.j  as  the  ferrite  material  js  allowed  to  become  isotropic.  Thus 
constructing  solutions  B^(a'kr)  and  H^(ajkr)  not  only  facilitated  the 
co.-nputationol  work  but  also  aided  in  physical  interpretation. 

Thus 

B^(0;x)  -  J^(x),  H^(0;x)  =  N^(x)  (3) 


where  x  =  kr. 

These  functions  are  given  by 


UB 

S,  (“ijx)  *  2  a  X 

1  ^  p-0  p 


P+1 


a^+2  -  -  (a  ap^^+ap)/(p+2)(p+)4) 


(9) 


and 


•1/2,  aj^  -  -t./3 


lL(c;x,  =  |{[Y  *  In  4))  B,(a;i) - k - 

^  ^  (l-abx  1*0“^ 


n+1 


*  Jo  *  * 


“  -  (a^^^-,+d^)/n(n+2)  -  [l/n  +  l/(n+2)]a 


n+1  n 


n 


(iO) 


..2 


d^  =  -  l/(l+a  ),  -  -  a/(l+a  ),  “  (r  -  Ip/  -  l/2)/2 


where 
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Figures  1  through  U  shov:  the  functions  B^(a;x),  H^(ajx)  and  their  der¬ 
ivatives  as  3  function  of  x  with  a  as  a  parai.jter  (-2  <  ®  <  10; 

<  X  <  10).  The  solution  to  equation  (3)  may  thus  be  written  in  the  form 

E.  =  AB-  Cajx)  +  BH,  (ajx)  (11) 


where  A  and  B  are  constants.  However,  E ,  must  be  bounded  at  the 

»  0 

origin  thus  B  “  0  The  boundary  conditions  (7)  then  gives 


h  (ajXo)  “  0 


iMiere  x^  ®  kb.  The  eigenvalues  of  (12)  will  be  denoted  by  (a)  where 
the  Correspond  to  the  eigenvalues  of  J^(x)  =  0.  Figure  5  shows 

the  functional  dependence  of  ^he  first  three  eigenvalues  p  .(u),  p.  ^(a) 

1,1  IfC 

and  p.  (a)  on  a. 

To  calculate  the  propagation  constant  at  given  co  and  given  ferrite 
parameters  we  use  relations  (U),  eliminate  P  to  obtain 


x^Ca) 


0 


x^(o) 


l+[a(l+x)Al^  l+a^(l+)(f/x^ 


Superposition  of  the  c\irve3  generated  by  (13)  on  the  family  of  oorves 
Pt  C^)>  give,  for  any  one  value  of  m,  the  two  sets  of  values  (a,x)  that 

aimiiltaneously  satisfy  the  required  conditions  (12)  and  (13).  Then,  using 
(4),  \  ^  may  calculate  the  propagation  constants,  i*e., 

pb  =  .  ax,  (lU) 

It  is  clear  that  the  curves  generated  by  (13)  are  symmetric  about 
a  =  0  and  the  cui’ves  p^^  ^(a)  are  not.  Thus  the  propagation  constants 
and  ^  calculated  from  their  points  of  intersection  differ.  We  observe 
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that  as  w  is  decreased  the  guide  is  'cut-off  in  one  direction  before  the 
other  and  the  ordinate  of  (13)  approaches  p,  (0),  Also,  since  the  slope  of 

•I* 

p  (a)  at  a  =  0  is  noirotive  and  the  slope  of  (13)  is  zero  at  a  =  0  -  we  see 
that  just  prior  to  'cut  ‘f  for  both  directions  of  propagation  there  are 
two  intersections  of  (13)  with  p.  (o)  for  a  positive  and  therefore  two 
distinct  modes  for  bhe  same  direction  of  propagation  become  possible.  These 
then  coalesce  as  ro  is  further  reduced  and  complete  cut-off  is  finally- 
reached.  The  directions  of  propagation  arc  not  identified  es  forw^^rd  und 
reverse  since,  for  a  positive,  |3b  may  be  positive  or  negative  depending  on 
the  sign  of  (I  X)/v.,  i.e.. 


^  [  (nO^  +  (i  x)/k  ^  o. 


US) 


V/e  can  noi:  treat  the  case  of  a  circumfei'entially  magnetized  ferrite  tube 
mounted  coaxially  vdthin  a  circular  cylindrical  g\iidc.  (We  again  restrict 
ourselves  to  considorin;  the  rotationally  symmetric  TE  modes).  Let  the  inner 
and  outer  radij,  of  -ihr  ferrite  tube  be  a  and  b,  respectively,  and  let  d 
be  the  guide  radius,  A  cross-section  of  the  guide  will  then  show  three  regions, 
namely,  an  inner  region  filled  by  some  isotropic  dielectric,  the  middle  region 
containing  the  ferrite,  and  an  outer  region  again  containing  an  isotropic 
dielectric,  those  regions  and  related  quantities  will  be  identified  using  the 
Roman  numerals  I,  II  and  III  for  the  inner,  middle  and  outer  regions. 

We  can  nov;  write  the  following  expressions  for  the  electric  field  com¬ 
ponents 

.J"  -  [i3^(a;kr)  +  CH^(a;kr)}  e"^^!!* 


(16) 


where 


’>  2 
k"  =  (a  ki>  fc 
o  o 


(17) 


and  where  ^  ^  constants.  The  dielectric  constant 

of  the  regions  not  Ctiotaining  ferrite  was  taken  to  be  ttiat  of  the  ferrite. 

This  m^,  in  fact,  be  desirable  in  practice  besides  sii-iplifying  the  calculations 
somewliat  • 


Application  of  the  boundary  conditions  and  eliminating  the  constants  yields 
the  transcendental  equation 


Bi(ajY‘)  8^(0;^) 

*  F3(aj  * 


Fi(a;  y<^)-F2(®; '2 
F^(a;  r^xh?^(a:  t^) 


(18) 


where 

F^CajTx)  =  [l-aTx+  tx  B|(a;Tx/B^(aj  Tx)j/Ci» 
F2(c;t>^  =  ly  JqC'V)/J]_(t  ) 

F^(a;  rx)  -  -aT>:  +  TxHj^(aj Tx;/Uj^/a;Tx)J/A‘ 

F^(a;'^^)  “Ty 

and 

ZQ(Ty)  =  JQ(Ty)  Nj^(y)  -  J^(y)  N^(Ty) 

Z^(Ty)  =  Jj^(Ty)  N^(y)  -  J^(y)  N^(Ty) 

where  t  may  te’-^e  on  the  subscripts  1  and  2,  and 

"  a/d,  Tg  «  b/d 

and 

X  =  kd,  y  =  k^d 


Si 


b( 

f< 


tx 


an 


iilSO 


'  .^y  ~2_  ^ ^  '-uC, 


We  aave  also 


2  - 
_  f .  w  .  r 


y'-  =  X  {x"-^La(ln)j"(l^A')}/v  4i'  (19) 

iolutj.on^  wf  o  <Q  tn^nvjGndental  equabio'i  v’all  yield  the  eir:ea\^.lues  as  a 
function  of  a  -  the  intersection  of  these  curves  with  (19)  will  then^  as 
before,  allow  computation  oi  the  propagation  constants  and  the  dif¬ 
ferential  pnase  shift  of  this  structure. 

a  sample,  and  cecause  these  results  v/ere  used  in  o.  iculations  to  be 
presented  below,  the  fom  that  the  fielc  components  take  for  the  ferrite 
tuuS  in  contact.  vcLuii  the  guide  wall  i?  given. 

For  0<r<a  \je  have 

%  =  J;,(lT/b)/j,(Ty) 

b  HJ  =  y  J^(yr/b)/j^(Ty)  (2o) 

b  hJ  =  jax  (1  +  X)  J^(yr/b)/x  J^(Ty) 


and  for  a  <  r  <  b  =  d 


^0  ~  F(ajxr/b) 


_  B^( a;xr/b  )H^( a;  x)-B^(  a;x)H^ ( a;  xr/b ) 
Bj^(ajTx)H^(ajx)-B^{a;x]H^(ajTx) 


-jcciA^b  -  {[b/r  -  ax]  F(a^xr/b)  +  x  ?' (a^xr/b) J /^' 

^  i[ ax(?i. -^xy^^  “  xb/(l  +X)r]  F(a;xr/b)  +  x  /(l  ^  x)» 


F'Cajxr/b)  , 


where  the  primes  . ndicate  differentiation  w.r.t.  xr/b. 


Henults  and  discussion 

As  was  to  be  expected  for  tiie  case  of  the  completely  filled  guide  the 
curves  (Fig.  5)  displaying  the  behavior  of  the  eigenvalues  as  a  function  of 
a  are  not  synmetric  about  the  ordinate  since  the  nagnetic  intensities  in 
the  regions  vine  re  circular  polarization  occurs  differ,  being  somewhat  stronger 
nearer  the  center  of  the  guide.  Therefore,  even  for  the  complete!;/  filled 
guide  a  differential  phase  shift  is  obtained.  This  is  in  contrast  to  the  case 
of  the  completely  filled  rectangular  guide. 

Figures  6  and  7  show  the  results  of  calculations  performed  on  the  trans¬ 
cendental  equation  (l6).  That  is,  they  show  the  behavior  of  the  eigenvalues 
as  a  function  of  a  for  the  two  limit  cases,  the  ferrite  tube  collapsed  into 
a  central  rod  and  the  tube  in  contact  with  the  guide  wall.  These  calculations 
were  made  taking  the  values  X  =  0  and  x  =  .2,  i.e.,  the  magnetizing  field 
H  was  taken  to  be  zero  and  rli/(o  was  given  a  representative  value  of  one 
fifth.  The  lowest  order  eigenvalue  was  calculated  for  various  values  of  t 
where  t  =  a/d  -  thus  t  indicates  the  fraction  of  the  guide  cross-section 
filled  with  ferrite!  In  the  limit  as  t  approaches  zero  or  unity  these 
cuives  approach  those  for  the  guide  filled  completely  by  an  isot-opic  dielectric 
or  the  ferrite.  The  dj electric  constant  for  both  the  ferrite  and  the  dielectric 
regions  was  tal-sn  to  be  nine. 

One  would  expect  that  in  the  partially  filled  guide  the  differential  phase 
shift  would  increase  over  that  of  the  comple+'^ly  filled  guide  as,  in  this 
case,  only  one  region  where  circular  polarization  occurs  is  contained  within 
the  ferrite  region.  Further,  one  expects  that  as  t  increases  the  dif- 
ferentialy  phase  shift  increases  until  t  is  such  that  about  half  the  guide 
section  is  filled  whereafter  it  would  decrease  until  it  matches  the  value 


for  the  com-letol,;  rolled  fr^ide.  Thus,  one  cxo-cts  the  slopf  of  the 
sigenvalue  cnrv'S  tc  be  the  largest  (negativel-)  uhen  t  is  about  otie 
half.  This  is  indeed  true  for  a  small  but  as  a  increases  negativclv 
this  no  lonser  holds  and  the  curves  take  a  dounturn  such  that  in  this  range 
the  differentia^  phase  shift  is  less  than  xt  would  be  for  the  con’'lately 
■guide,  j.  :n  u’e  have  a  ferrite  tube  iti  contu.ct  v.’jt!;-  t  he  gi^iid.i  ’.'cll 
this  effect  is  oven  more  pronounced  (Fig,  7)  and  ue  see  that  for  oven  smuII 
u  the  slopes  of  tli'O  eigonralus  c'lrves  are  s-iurl  to  or  less  than  that  for 
the  complete uy  i^llcci  {Tiiide,  Of  course,  for  the-  latter  case  v;here  the  guide 
contains  little  lerrite  ( 1.)  we  exooct  t-ivc  eigenvalues  carves  tc 
become  almost  sycnnetric  about  the  ordinate  since  the  differential  phase  shift 
will  become  vanisninglgv  small. 

An  explanation  that  would  account  for  the  dox^nturn  in  the  eigenvalue 
curve  for  a  large  and  negative  is  that  the  field  muse  bo  exnell-^d  from 
the  ferrite  in  this  region.  To  determine  if  in  fact  this  occurs  tne  com- 
ronents  of  the  electric  and  magnetic  field  v/ere  computed  for  this  case 
using  (20)  and  (>1)  for  .2  and  x^Co)  =  6,  20,  30  and  hO  in  (13). 

These  values  of  "^g(O)  inserted  into  (13)  assured  that  the  intersection 
points  of  (13)  ■  ith  tne  eigenvalue  curve  for  t  =  .2  fell  at  suitcble  points, 
that  is,  boforo.  on,  and  after  the  downturn.  The  results  cf  these  calcula¬ 
tions,  presented  in  Figure  3,  clearly  show  that  the  fields  are  excelled  as 
®  increases  n.u'atjvcly  —  until  at  Xq(0)  eoual  to  alxjut  ;0  the  greater  oart 
of  the  ener.gx'  is  transported  for  this  direction  of  prou-agotion  in  the  isotro¬ 
pic  region.  Thui,  for  this  case  the  magnetic  fields  are  relatively  v/eax 
throughout  the  region  of  circular  polarization  IjnLng  within  the  ferrite  and 
conseviuently  tr.u  ar.j sotropio  rpj.racter  of  tne  medium  htis  a  smaller  effect  on 


the  nronagaticn  constant  for 


of  propagation. 


this  direction 
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The  region  of  a  large  and  negative  holds,  ia  fact,  but  j^odest  intarost 
since  operation  in  this  region  would  yield  decreasing  values  of  differential 
phase  shift.  Thus  the  character  the  eigenvalue  curves  for  a  above  the 
downtown  is  of  nrjnar;/'  interest. 

Figure  9  presents  a  typical  curve  of  the  differential  phase  shift  as  a 
function  of  t.  Again,  the  explanation  of  the  unexpected  result  that  the 
maximum  differential  phase  shift  occurs  for  the  completely  filled  guide  lies 
in  tte  fact  that  the  fields  show  a  preference  for  the  isotronic  over  the 
anisotropic  region,  Vfe  note  that  when  the  guide  is  about  tv;o  thirds  filled 
with  ferrite  ('f  =  .36)  the  effect  introduced  by  the  region  of  circular 
polarization  near  the  guide  is  exactly  counterbalanced  by  the  effect  felt 
due  to  the  ferrite  approaching  the  region  of  circular  polarisation  near  the 
guide  axis. 

The  final  ciilculations  presented  here  were  made  to  observe  the  changes 
in  the  differential  phase  shift  for  a  thin  tube  of  ferrite  into  the  guide 
for  various  values  of  the  tube  mean  radius.  The  computations  vTQre  made  for 
a  ferrite  tube  of  thickness  one  tenth  the  guide  radius,  the  parameters  X 
and  X  again  being  0  and  .2,  respectively,  while  the  remaining  guide 
rarameters  were  selected  such  that  ^^(0)  =  6.  The  calculations  of  the  dif¬ 
ferential  phase  shift  variation  with  mean  tube  radius  are  shewn  in  Figure  10. 
This  graph  clearly  demonstrates  that  the  differential  phase  shift  Increases 
until  the  ferrite  tube  contains  a  region  of  circular  polarization  tnen 
decreases,  passes  through  zero,  and  decreases  to  a  negative  maximum  when  the 
tube  encloses  the  region  circularly  polarized  in  the  opposite  sense.  The 
fact  that  the  maximum  differential  phase  shift  is  tbe  greater  for  true  region 
of  polarisation  nearer  the  guide  axis  confirms  that  the  ma  - ?tic  i ntensity  is 
the  greater  there. 
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Even  though  the  results  for  the  thin  tube  may  be  arrived  at  through 
perturbation  theory  the  exact  solution  allows  the  extension  of  the  cal^ 
culationto  thick  tubes  and  to  the  investigation  of  the  effect  of  hf-hcr 
order  modes  which  may  yxeld  greatly  increased  differential  phase  shift. 


FIG. 
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